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(57) Abstract: The invention relates to a method of estimating a rainfall rate in the case of rain which only generates negligible 
attenuation using a bipolar radar. The inventive method comprises the following steps: the bipolar radar is used to measure the 
differential phase $ DP and the apparent reflectivity Zc along at least one of the polarisations H or V over a given segment [r 0 , rj of a 
ray path in relation to said radar; an estimate of value N 0 *, which is representative of the dimensional distribution of the rain drops, is 
determined using the difference in differential phase between r 0 and r! and using an integral of a function of the apparent reflectivity 
Zc along the length of segment [ro, rj; the precipitation rate value is derived at a point using N G * and the apparent reflectivity at said 
point. 



S (57) Abrege : L'invention concerne un precede' d' estimation d'un taux de precipitation dans le cas d'une pluie ne generant qu'une 
attenuation negligeable, au moyen d'un radar bipolaire, consume par les Stapes suivantes: on mesure, au moyen dudit radar bipolaire, 

2 sur un intervalle [ro, rj donne de rayon de trajet par rapport audit radar, la phase differentielle * DP et la reflectivite apparente Ze selon 
au moins Tune des polarisations H ou V; on determine une estimation de la valeur N 0 * representative de la distribution dimensionnelle 

^ des gouttes de pluie, a partir de la difference de la phase differentielle entre r 0 et ri et a partir d'une integrate d'une fonction de la 
r6flectivit6 apparente Ze le long de V intervalle [r 0 , T\]\ on d^duit la valeur du taux de precipitation en un point a partir de N D * et de la 

^ reflectivite" apparente en ce point. 



WO 03/007016 PCT/FR02/02447 

1 

« Technique pour Testimation de pluie avec radar lorsque I'attenuation 
est negligeable » 

La presente invention est relative aux techniques m6teorologiques 
pour Pestimation d'un taux precipitant au moyen d'un radar. 
5 ^estimation du taux de pluie a partir de mesures de radar est 

g6n6ralement confrontee au probleme de ('attenuation de I'onde radar, et a 
celui de la variability naturelle de la pluie. 

Le concept d'un radar coherent & diversity de polarisation associe a 
un algorithme appele ZPHI a ete decrit dans WO 99/38 028, comme solution 
10 apportee a ces deux obstacles dans des conditions op6rationnelles. 

ZPHI est un algorithme profileur qui, & la base, utilise en entree un 
profil de reflectivite mesur6e Za et une contrainte donnee par une difference 
de phase differentielle O dp entre deux points n et r 2 d'une ligne de visee. 

A partir de ces mesures, on determine I'attenuation specifique A et un 
15 parametre connu sous le nom de N 0 \ parametre quantitatif de la distribution 
des tallies de gouttes. 

Le taux de pluie R a estimer est obtenu comme une fonction de ces 
deux parametres, pour les deux raisons suivantes : 

- I'attenuation specifique A n'est pas soumise aux effets 
20 d'attenuation. Son utilisation dans Testimation de R permet ainsi de 

regler le probleme de I'attenuation. 

- le parametre N 0 * suffit pour dScrire la variabilite naturelle de la 
pluie. 

L'algorithme ZPHI s'applique aussi bien en bandes X et C, bandes de 
25 frequence sensibles a I'attenuation. 

On vise ici maintenant un proc§d6 6tendu a partir du procede ZPHI, 
lorsque I'attenuation n'est pas importante. 

Pour ce cas particulier, il est possible de prescrire un formalisme base 
sur des elements de base, mais plus adapte et beaucoup plus direct. 
30 L'algorithme ZPHI peut 6galement etre applique a la bande S, pour 

laquelle I'attenuation n'est pas importante, voire negligeable. En effet, son 
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interet subsiste dans la mesure oil il reproduit « en direct » la variabilite 
naturelle de la pluie. Le recours au parametre A n'est alors pas forcement 
necessaire. 

L'invention a done pour but de proposer un procede simple, fiable et 
5 efficace devaluation du taux precipitant dans le cas d'attenuation faible ou 
negligeable. 

L'invention propose pour cela un procede d'estimation d'un taux de 
precipitation dans le cas d'une pluie ne generant qu'une attenuation 
negligeable, au moyen d'un radar bipolaire, constitue par les etapes 
10 suivantes : 

- on mesure au moyen dudit radar bipolaire, sur un intervalle [r 0 , n] 
donne de rayon de trajet par rapport audit radar, la phase differentielle <D D p et 
la reflectivite apparente Z e selon au moins I'une des polarisations H ou V ; 

- on determine une estimation de la valeur N 0 representative de la 
15 distribution dimensionnelle des gouttes de pluie, a partir de la difference de la 

phase differentielle entre r 0 et n et a partir d'une integrate d'une fonction de 
la reflectivite apparente Z e le long de I'intervalle [r 0 , r] ; 

- on deduit la valeur du taux de precipitation representant en tout point 
a partir de No* et de la reflectivite apparente en ce point. 

20 On propose egalement selon l'invention un dispositif pour I'estimation 

d'un taux precipitant comportant un radar bipolaire, ainsi que des moyens de 
traitement, ledit radar comportant des moyens pour la mesure de la phase 
differentielle et de la reflectivite selon au moins une des polarisations H ou V, 
caracterise en ce que les moyens de traitement mettent en ceuvre les 

2 5 differentes etapes de traitement du procede selon le paragraphe precedent. 

D'autres caracteristiques, buts et avantages de l'invention apparaitront 
a la lecture de la description detaillee qui va suivre. 

Comme pour le procede ZPHI, le procede decrit maintenant requiert 
un radar a diversite de polarisation coherent. Les donnees d'entree sont le 

30 profil de reflectivite Z H (ou Z v ), ainsi que la mesure de la phase differentielle 
entre la voie H et la voie V. 



WO 03/007016 



PCT/FR02/02447 



3 



On vise ici a determiner le parametre N 0 * (ainsi que le taux de pluie) 
dans les cas ou I'attenuation est negligeable, par exemple cas de la pluie 
mesuree par la bande S, cas de pluies suffisamment faibles pour les bandes 
C et X, ou encore cas de la glace. 
5 Comme dans le procede ZPHI, le present procede integre deux types 

de donnees : la reflectivite Z H (ou Z v> ou plus generalement Z e ) et la phase 
differentielle cj> DP . 

Dans la mesure ou I'attenuation est negligeable, la reflectivite 
mesuree peut etre directement utilisee dans I'estimation du taux de pluie R, 
10 moyennant une estimation du parametre N 0 * qui est precisement I'objet du 
present procede. 

Le procede repose sur la relation « universelle » qui lie la reflectivite 
equivalente Ze (mm 6 m' 3 ) et le taux de phase differentielle K D p(°km" 1 ) : 



15 



20 



K 



DP 



= a 



L*o J 



(1) 



Ou a et b sont des coefficients specifies par le modele de diffusion, qui 
dependent du type de precipitations (pluie ou type de cristal de glace) et qui 
sont legerement fonction de la temperature. La determination du type de 
precipitations peut etre operee par une methode de classification du type de 
celle decrite par Straka, Zrnie et Ryzhkov (2000). 

Plus generalement, le modele de diffusion definit la relation : 



K 



£L = F 



(2) 



En integrant (1) ou (2) entre les deux bornes r n et r 2 du segment 
d'integation (^<t 2 ), on obtient : 



25 



0 DP (r 2 )-0 D/) (r 1 ) = [<] 1 - t j; i 2 z/* 



(3) 
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Soit: 



N n = 



a£z e b ds 



1-6 



(4) 



II ne reste plus alors qu'a transporter cette estimee de No* dans la 
relation « universelle » liant R a Z e : 



10 



15 



20 



25 



R = cM-%" 



(5) 



pour obtenlr I'estimee de R. 

Dans le cas general ou la relation K D p-Ze s'ecrit sous la forme (2), (3) 



devient 



ds 



(6) 



On n'a alors pas de solution analytique pour N 0 \ mais une solution 
numerique peut etre trouvee par methode numerique iterative en utilisant 
comme valeur de depart pour N 0 *(« first guess ») la solution (4). 

Dans le cas de la glace, qui n'attenue pratiquement pas I'onde radar 
quelle que soit la frequence utilisee (bandes X, C et S), les limites 
d'application ne sont pas fixees. 

Dans le cas de la pluie, I'application du procede est contrainte par une 
attenuation negligeable entre les points ^ et r 2 (r 2 >ri). Plus specifiquement, 
les conditions d'application qui conduisent a une erreur de 3 dB sur N 0 sont 
(pour T = 10°C, et une loi de distribution de failles de gouttes de type gamma 
avec un parametre de forme egal a 2) : 

Pour la bande X : [<D D p(r 2 ) - 0 D p(ri)] = <, 4° 

Pour la bande C : [(Dopte) - <J>Dp(ri)] = < 1 0° 

Pour la bande S : [<J> D p(r 2 ) - Owfri)] = ^ 64 ° 
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Les domaines d'application principaux de I'invention sont : 

les memes domaines d'application que dans le cas de 
I'algorithme ZPHI, a savoir I'estimation de la pluie sur les 
bassins versants pour la surveillance des crues et la gestion 
5 des ressources en eau. Cette application est valable pour tout 

type de pluie pour la bande S, ou aux pluies suffisamment 
faibles pour les bandes C et X. 
I'estimation des precipitations en phase glace. 
L'invention presente egalement de nombreuses autres applications en 
10 meteorologie. 

[1] Differential propagation phase shift and rainfall rate estimation, 
Sachidananda M. et Zrnie D.S., Radio Science, 21-2, p.235-247 
(1986). 

[2] Polarimetric method for ice water content determination, Ryshkov 
15 A.V., Zrnie D.S. et B.A. Gordon. J. Appl. Meteor., p1 25-1 34 (1998). 

[3] The rain profiling algorithm applied to polarimetric weather radar, 
Testud J.E. Le Bouar, E. Obligis et M. Ali-Mehenni, J. Atmos. Oceanic 
Technol., 17, p332-356 (2000). 
[4] Bulk hydrometeor classification and quantification using polarimetric 
20 radar data: Synthesis of relations, Straka J.M., D.S. Zrnie, A.V. 

Ryzhkov, J. Appl. Meteor., 39, p1 341-1372 (2000). 
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REVEND1CATIONS 



1 . Procede d'estimation d'un taux de precipitation dans le cas d'une 
5 pluie ne generant qu'une attenuation negligeable, au moyen d'un radar 

bipolaire, constitue par les etapes suivantes : 

- on mesure, au moyen dudit radar bipolaire, sur un intervalle [r 0 , ri] 
donne de rayon de trajet par rapport audit radar, la phase differentielle <D D p et 
la reflectivite apparente Ze selon au moins Tune des polarisations H ou V ; 

10 - on determine une estimation de la valeur N 0 * representative de la 

distribution dimensionnelle des gouttes de pluie, a partir de la difference de la 
phase differentielle entre r 0 et n et a partir d'une integrale d'une fonction de 
la reflectivite apparente Z e le long de I'intervalle [r 0 , n] ; 

- on deduit la valeur du taux de precipitation en un point a partir de No* 
15 et de la reflectivite apparente en ce point. 

2. Procede selon la revendication 1, characterise en ce qu'on deduit N 0 * 
directement a partir de la difference de phase differentielle entre r 0 et n et a 
partir de I'integrale de la reflectivite apparente Z e elevee a un exposant 
choisi. 

20 3 .Procede selon la revendication 2, caracterise en ce que I'exposant 

choisi est un exposant b, qui verifie : 




ou K D p est le taux de variation de la phase differentielle le long du rayon, a et 
b etant specifies par le type de precipitation considere et par la temperature. 
25 4. Proced6 selon la revendication 3, caracterise en ce que le type de 

precipitation est un type correspondant parmi le type pluie et le type cristal de 
glace. 
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5. Precede selon la revendication 1 , caracterise en ce que la fonction 
de la reflectivite apparente qui est integree est dependante de N 0 \ et en ce 
que le procede inclut le fait de calculer a plusieurs iterations cette fonction, 
avec une valeur No* determinee a ('iteration precedente. 

6. Procede selon la revendication 5, caracterise en ce que Ton calcule 
N 0 * a partir de la relation : 



2 ( Z } 



10 



ou Odp est la phase differentielle et oil F est une fonction verifiant la 
relation : 



K 



£L = F 



ou Kdp est le taux de variation de la phase differentielle le long du rayon. 

7. Procede selon la revendication 5 ou la revendication 6, caracterise 
15 en ce que I'on utilise, a une premiere iteration, la valeur de N 0 * donnee par : 



4, Z > 



1-6 



8. Procede selon I'une quelconque des revendications precedentes, 
caracterise en ce que Ton determine le taux de pluie R a partir de N 0 * et de la 
20 reflectivite apparente Z e par la relation : 
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9. Dispositif pour ('estimation d'un taux de precipitation comportant un 
radar bipolaire, ainsi que des moyens de traitement, ledit radar comportant 
des moyens pour la mesure de la phase differentielle et de la reflectivite 
selon au moins une des polarisations H ou V, caracterise en ce que les 
5 moyens de traitement mettent en oeuvre les differentes etapes de traitement 
du procede selon Tune quelconque des revendications precedentes. 
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